Zero and low temperature behavior of the two-dimensional ±J Ising spin glass 
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Scaling arguments and precise simulations are used to study the square lattice ±J Ising spin 
glass, a prototypical model for glassy systems. Droplet theory explains, and our numerical results 
show, entropically-stabilized long range spin glass order at zero temperature, which resembles the 
energetic stabilization of long range order in higher dimensional models at finite temperature. At 
low temperature, a temperature-dependent crossover length scale is used to predict the power-law 
dependence on temperature of the heat capacity and clarify the importance of disorder distributions. 
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Glassy systems, characterized by extremely slow re- 
laxation and resultant complex hysteresis and memory 
effects, are difficult to study because their dynamics en- 
compass a great range of time scales [T] . Glassy materi- 
als include those without intrinsic disorder, such as sil- 
ica glass, and those where quenched disorder influences 
the active degrees of freedom. An example of the lat- 
ter is the Edwards- Anderson spin glass model [Ij, which 
includes the disorder and frustration necessary to cap- 
ture many of the complex behaviors seen in disordered 
magnetic materials. Many aspects of this prototypical 
spin glass model remain poorly understood. The droplet 
and replica-symmetry-breaking pictures provide distinct 
views of spin glasses [3]. Analytical results are few, so 
numerical approaches are invaluable for both testing and 
motivating new ideas. But numerics can also be exceed- 
ingly difficult in the general case: computing spin glass 
ground states is believed to require exponential time to 
solve exactly, as it is an NP-hard problem [U [5] . 

A fortunate special case that avoids this computa- 
tional intractability is the two-dimensional Ising spin 
glass (2DISG). The Hamiltonian is H ~ JijSiSj. 
The nearest-neighbor couplings — {Jij} are indepen- 
dent random variables coupling Ising spins = ±1 at 
sites j on a square lattice with sites. The random signs 
of the Jij lead to competing interactions that can not be 
simultaneously satisfied. Here we study the model with 
the ±J distribution, where each bond value is Jij — ±1 
with equal probability. Highly developed numerical al- 
gorithms ^ EHH] can efficiently circumvent the complex- 
ity due to disorder and frustration: ground states and 
finite-temperature partition functions of the 2DISG may 
be computed in time polynomial in L. These algorithms 
have given us more insight into model glassy systems. 

In this Letter, we develop the droplet scaling theory as 
applied to the zero and low temperature properties of the 
± J model. We obtain reliable numerical results for very 
low temperatures T and large L (well beyond any Monte 
Carlo simulation results), including the probability dis- 
tributions of spin correlations. We show that the ±J 
model at T = in 2D has equilibrium correlations simi- 



lar to those of a higher-dimensional spin glass at T > 
in its spin glass ordered phase. Thus the 2DISG provides 
a computationally-tractable model with spin-glass order 
that closely mimics spin glass order in higher dimensions, 
where large samples can not be studied numerically. 

The ground state of a ± J model is highly degenerate 
with an extensive entropy [6l [TOl [11] . For many years it 
was assumed that the ground state is critical, with power- 
law spin-spin correlations as a function of distance |111 - 
[13]. Jorg, et al. [TS] then presented evidence that instead 
the ground state has true long-range spin glass order. Re- 
cent results [16] about rigid clusters of spins with fixed 
relative orientations are also suggestive of long range cor- 
relations. Here we add to that evidence, and present the 
corresponding droplet theory, which differs in some im- 
portant respects [T7] from the scenarios proposed in Refs. 
[I5] [TBI [18] . We show how the T = properties predict a 
new temperature-dependent crossover length scale ix{T) 
where the probability distribution of the droplet free en- 
ergies crosses over from discrete to continuous. We use 
this to predict the low-T scaling of the specific heat. This 
prediction is shown to agree well with our numerical re- 
sults, resolving the long-debated question of the specific 
heat behavior in this model. These results provide new 
insight into the role of entropy and the importance of the 
choice of disorder distribution in spin glasses. 

Computational technique We start by computing 
the partition function Zp{J^) in an L-hy-L sample with 
periodic boundary conditions, adopting previously pub- 
lished techniques [HI [7], but using arbitrary precision 
arithmetic [5]. The efhcient computation of Z relies on 
a hierarchical decomposition of the sample 7J: start- 
ing from the smallest pieces of the spatial decomposi- 
tion, single plaquettes, neighboring pieces are joined to- 
gether to recursively compute the partition function. The 
same four Pfaffians JSj used to find Zp also give Zap, 
the partition function for antiperiodic boundary condi- 
tions, where the horizontal bonds along a vertical col- 
umn have negated Jij. Note that these computations 
are exact to within numerical precision: there is no ques- 
tion of convergence as there is with Markov chain Monte 
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Carlo methods and we have verified the numerical sta- 
bility of our procedure [9]- The two partition functions 
differ due to a domain wall induced by the change in 
boundary conditions. We use finite differences over T 
of F = —kTlii{Z) to compute equilibrium averages for 
energies E and entropies S. The sample-dependent dif- 
ferences dX{J) = Xp - Xap, with X = F, E, and S, 
give relative domain wall (free) energies and entropies. 
Error bars in our plots indicate la statistical errors. 

We have extended this technique to compute arbitrary 
spin-spin correlation functions [H]. These correlations 
can be computed at sufficiently low T, T — 0.02 for 
L < 256, and with sufficient precision, 4096 bits, such 
that the contribution of excited states can be clearly sep- 
arated from the ground state contribution to correlations. 
If spins have fixed relative orientation with probability 
greater than 1 — e"^/"^, we take the spins to be rigidly 
correlated at T = 0. Such spins have fluctuations that 
are clearly smaller than non-rigid spins by a factor of at 
least « 10^^, exceeding multiplicity (entropic) factors for 
excited states. We have confirmed that the assignment of 
rigid correlations is unchanged if T is lowered or the pre- 
cision increased, using over 200 samples for L — 128. An 
example of T = rigid and non-rigid nearest-neighbor 
spin correlations (siSj)o, and the domain wall due to a 
change in boundary conditions are shown in Fig. [l] 




FIG. 1: [color online] An L — 128 sample, (a) Nearest 
neighbor correlations: Black lines are bonds that are rigidly 
correlated at T = 0. (b) Relative domain wall. Black lines 
are bonds whose correlations are not affected by the change 
P — >■ AP in the horizontal boundary conditions. The shaded 
bonds indicate the change in the correlations on a logarithmic 
scale, with light colors for bonds with \S{siSj)o\ ~ 1, and the 
darkest colors indicating \5{siSj)o\ ~ e"^". This "diffuse" 
domain wall appears to have a multifractal geometry [20j . 

Domain walls, droplets, and rigid clusters Con- 
sider a given ground state configuration of a particular 
realization J' of the ±J model. A relative domain wall 
loop is the boundary of a simply connected region where 
all the spins are flipped relative to the given ground state 
configuration. The droplet excitation [21] at site i and 
length scale £ is the set of all lowest-energy domain wall 
loops of linear size of order £ enclosing site i and of order 
£^ other sites. We can ask about the probability distribu- 
tions of the energy and of the entropy of this droplet ex- 



citation. Here we do not directly measure these droplets, 
but computations for the very similar system-size rela- 
tive domain walls due to changing between P and AP 
boundary conditions give insight into the effect of these 
droplets on the spin-spin correlations. 

We flrst note previous work (e.g., which shows 

that system-size domain walls have a probability distri- 
bution for SE that is independent of L, as L — >■ oo (given 
a parity for L). For large even L, we flnd that 77. 5 ±1.0% 
of the samples do have zero energy domain walls. Thus, 
assuming the droplets have a similar probability distri- 
bution, the event that there is no zero energy droplet of 
scale i surrounding a site i has probability strictly be- 
tween zero and one. This is also true for smaller scales 
£/2, £/A, etc., with the correlations between the events 
at each scale presumably falling off rapidly across scales. 
As a result jl4j , a given site i has no zero energy droplets 
at any scale less than £ with a probability p{£) that de- 
creases with increasing £ as a power of £. 

In a finite sample at T = 0, we choose to identify 
the "backbone" of the ground states as the largest set of 
rigidly correlated spins |16j . Spins on the backbone can 
not be flipped at zero energy by any droplets without 
flipping the entire backbone. By the above argument, 
the number of such spins scales with a power of L less 
than two: the backbone is a fractal. Hartmann's results 
imply a backbone fractal dimension di, = 1.78(2) [T3]; our 
results for rigidly correlated clusters are consistent with 
this (20j . As seen in previous work [TB] , we sometimes see 
backbone spins that are not connected by a path of rigid 
bonds; the backbone can have "gaps". Thus the back- 
bone of the sample in Fig. 1(a) is all spins on the largest 
cluster connected by rigid bonds plus possibly some other 
spins or clusters whose rigidity is not detected by nearest- 
neighbor correlations. The backbone on this sample per- 
colates; the fraction of large samples where the backbone 
percolates is a number between zero and one [20] . 

Entropic enhancement of correlations So far, we 
have discussed the rigid correlations that arise from lo- 
cations where there are no zero energy droplets. What 
about zero energy droplets, where spin clusters flip in the 
ground state? Saul and Kardar [6] found that zero-energy 
relative domains walls have a typical entropy 5S ~ L^^ , 
with 6s = 0.50. Our data from L = 32 to L = 256 is 
very well fit by a power law, with ^5 — 0.50 ± 0.01, (in 
contrast with Ref. [23]). It is intriguing that this expo- 
nent is consistent with the simple rational number 1/2. 
The probability distribution of 6S is continuous through 
zero, so \5S\ is less than or of order one with a probability 
~ L~^^ . We take the probability distribution for 6S for 
the zero-energy droplets to have the same scaling. 

This leads to an essential new aspect of the droplet pic- 
ture for the two-dimensional ± J Ising spin glass at T = 0: 
A droplet at a large £ typically either has a nonzero en- 
ergy so never flips or has a large entropy so it is flipped 
away from its usual orientation only with exponentially 
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small probability ~ exp (— It is only the entropically 
active droplets with \6S\ less than or of order one that 
flip at T = with a substantial probability. These active 
droplets are power-law rare in the length scale i: only 
a fraction ~ £~^^ of the sample's area is active at scale 
£. This scenario is very similar to the expected behavior 
within the droplet theory |21j for an Ising spin glass in 3 
or more dimensions at a nonzero temperature within the 
spin glass ordered phase, where only a fraction ^ of 
the droplets at scale i have SF less than order of order 
kgT so are thermally active, and the typical free energy 
of a droplet scales as 5F ~ with 9 > 0. 

This droplet picture allows us to apply many of the 
predictions of the droplet theory [21] , with the usual free 
energy exponent 9 replaced by the entropy exponent 9s — 
0.50. In particular, it predicts long range spin glass order, 
with the T — spin-glass correlation function for spins at 
spacing r, Go{r) = [{S(^sp)q], where the square brackets 
denote an average over samples, behaving at large r as 



Go(rO-Go((X3)^r- 
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There is long range order because the total probability, 
found by summing over all scales, of flipping a given spin 
is bounded away from unity, since the probability of the 
droplet being entropically active decays as £~^^ . The 
active droplets with scale £ > r are unlikely to flip only 
one of a pair of a spins separated by distance r, so such 
droplets do not substantially reduce Go{f). Thus Go(r) 
is larger than Go(oo) by an amount proportional to the 
probability ^ t-^^s duplets of scale r being active. 

The results of numerical calculations for Gt {r) for T <C 
L~^^, shown in Fig. 2(a), are consistent with Eq. (1). 
To verify this and to more precisely determine the long 
range order Go(c»), we measure Gt{t) for higher T. At 
nonzero temperature and finite L, the difference between 
GT(r) and Go(oo) should scale as L-^^ g{f/ L,TL^^), for 
a scaling function g. The scaling parameter TL^^ nat- 
urally arises from the competition between energy and 
entropy. From such scaling of the correlations for several 
values of r/i, we infer Go(oo) — 0.395 ±0.010. A specific 
example of a scaling collapse, with fixed 9s = 0.50 and 
variable Go(oo), is shown in Fig. [2jb) for r — {L/A,L/A). 
In fact this long range order with Go(oo) = 0.4 is appar- 
ent in Poulter and Blackman [I^ (their Fig. 6), although 
that is not how they interpreted those data. A similar 
entropically-stabilized long range order has recently also 
been seen numerically in a diluted ±J spin glass in 3D 
where the ground state stiffness is zero [24] . 

Crossover to effectively continuous disorder Us- 
ing this droplet picture, we can also understand the 
crossover to the T > behavior. The competition be- 
tween entropy and energy [ITl |2T1 [25l |26] in spin glasses 
causes extreme sensitivity of the long-distance spin cor- 
relations to changes in T ("chaos"). We use a related 
argument at very low T to derive the low-T scaling of 
the specific heat in the ± J model. 
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FIG. 2: (a) Correlation functions at low temperatures for 
L = 32,64,128,256 plotted vs. r'^'' , using Os = 0.50. 
Droplet theory based on the entropy of the droplets pre- 
dicts a linear curve at small values of r~^® and large L, 
as is seen, (b) Sample scaling collapse for correlation func- 
tion Grif)- Using nonzero-temperature scaling, we expect 
i*s[GT(r) - Go(oo)] ^ g{r/L,TL^s). Adjusting Go(oo) for 
the best collapse, we find Go(c>o) = 0.395 ± 0.010. 



At r = 0, the droplet energies 5E are integer multiples 
of 4, so the distribution of droplet energies is discrete. 
However, the droplet entropies 5S{T = 0) ^ t^^ are con- 
tinuously distributed for large £. We define a crossover 
length £x{T) to be the scale at which the typical T5S{£x) 
becomes 0(1), causing 5F — SE — TSS to become con- 
tinuously distributed. Thus at low T 
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For £ <ti £x (T) and very low T only zero-energy droplets 
can be active. But at scales at and above £x{T) the 
droplets with 6E = 4 or more can, due to their large 
entropy, have SF near zero and thus also be thermally 
active (Ref. 15 gave a discussion of this crossover from 
discrete to continuous behavior, but not of £x(T).) 

Specific heat The low-T specific heat G in a 2DISG 
is governed by the thermally active droplets (those with 
excitation free energies 6F < 0{T)) that have nonzero 
energy |21l , with the smallest droplets dominating, since 
they have the highest density. If the disorder distri- 
bution is continuous, these smallest droplets have size 
0(1), energy of order T and density proportional to T, 
so give C ~ T at low T. For ± J disorder, on the other 
hand, the droplets with size £ < £x{T) have either zero 
energy SE or SE ~ 0(1) so do not contribute to the 
low-T specific heat. The smallest active droplets with 
SE are of size of order £^(T). Due to their 0(1) 
energy gap, these droplets each contribute ^ 1/T^ to G. 
They occur with density ^ T/t'^{T) so the specific heat 
C ~ r(2/es)-i at low T in an infinite sample (alterna- 
tively, the entropy density is £l?^'^ , which scales as G). 
For finite L, this power-law specific heat is cut off at the 
lowest temperatures when the size ix (T) of these active 
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droplets exceeds L. This produces the low-T finite-size 
scaHng form C « r(^/^^^~^c(TL^^), where c(x) is a scal- 
ing function that goes to a constant for large argument 
(where 1 ^ tx{T) <C L). Our specific heat data are 
shown in this scaling form in Fig. [3j Our data for in- 
termediate temperatures, T = 0.35, are consistent with 
Monte Carlo results [37] , which saw an effective exponent 
a ^ —4.2, with C ~ T~". However, the effective expo- 
nent crosses over to a lower magnitude consistent with 
a = 1 — TT- — —3.0 at low T and the largest L. Contri- 

tis " 

butions from the largest active droplets give a subleading 
correction, which from standard hyperscaling is a' = —2v 
(though see Ref. [TT]). that is too weak to detect yet. Our 
result, a = 1 — gives the leading contribution to C at 
low T for the 2D ±J model in the infinite volume limit. 
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FIG. 3: [color online] Scaling plot for the low-T specific heat 
C. For T < L^^ , there are no active droplets, so C is exponen- 
tially suppressed. The low-T scaling behavior in large-enough 
systems is seen when CT^-^/^^ jg g^t at TL®s > 1. The value 
|28| of 6s = 0.50 is independently set by the scaling of the 
entropy of domain walls. 

Discussion The ±J spin glass model in 2D can be 
used to efficiently study the thermodynamics of a glassy 
model. Large scale simulations verify conceptually novel 
scaling arguments for this model: these concepts have 
broad implications. Strictly correlated spins at T = 
form fractal clusters that can be explained using a droplet 
theory [T3], leading to backbone correlations that decay 
as a power law. Yet the spin-spin correlations at T = 
have long range order. Though large zero-energy ex- 
citations are abundant, the scale-dependent entropy of 
these excitations prevents them from destabilizing the 
long range order. This entropic "stiffness" at T = is 
very similar to the (free-)energetic stiffness in d > 3 mod- 
els that gives a stable spin glass with long range order 
at nonzero T, making the 2D ±J model a good candi- 
date for exploring the spin glass phase in general. The 
difference between discrete and continuous disorder is a 
general problem in statistical mechanics; here we show 
how this difference can have striking effects on the scal- 
ing through a temperature dependent crossover length. 
When the temperature is nonzero, application of this 
crossover length ix{T) leads to a new prediction for the 



low temperature specific heat of the ±J model in 2D, 
which we have verified numerically. 
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